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Abstract 

Metaheuristic and self-organizing criticality (SOC) could contribute to robust 
computation under perturbed environments. Implementing a logic gate in a computing 
system in a critical state is one of the intriguing ways to study the role of metaheuristics 
and SOCs. Here, we study the behavior of cellular automaton, game of life (GL), in 
asynchronous updating and implement probabilistic logic gates by using asynchronous 
GL. We find that asynchronous GL shows a phase transition, that the density of the state 
of 1 decays with the power law at the critical point, and that systems at the critical point 
have the most computability in asynchronous GL. We implement AND and OR gates in 
asynchronous GL with criticality, which shows good performance. Since tuning 
perturbations play an essential role in operating logic gates, our study reveals the 


interference between manipulation and perturbation in probabilistic logic gates. 


1. Introduction 


As we confront a large amount of data, it is difficult to obtain a complete dataset, resulting in an 
incomplete and ill-defined problem. Therefore, metaheuristics are investigated more intensely than 
optimization (Erskine & Hermann, 2015). The first direction of metaheuristics is bioinspired 
computing based on swarm intelligence, such as the krill herd algorithm (Wang et al., 2014), monarch 
butterfly optimization (Chakrabarty et al., 2014), ant colony algorithm (Dorigo & Stiitzle, 2004), and 
others (Wang, Deb & Coelho, 2018). These algorithms can search for quasi-optimistic solutions by 
balancing the global goal of a group with internal perturbations. The second direction is the 
development of biological computing, in which computing is implemented by living biological and/or 
chemical material, such as BZ pattern computing (Toth et al. 2009), physarum computing (Nakagaki, 
Yamada & Toth, 2000; Takamatsu et al., 2001; Tero et al., 2010; Tsuda, Aono & Gunji, 2004; Gunji 
et al., 2008; 2011; Adamatzky, 2010; Aono et al., 2015; Zhu et al., 2018), fungal computing 
(Adamatzky, 2018; Adamatzky et al., 2020) and soldier crab computing (Gunji, Nishiyama & 
Adamatzky, 2011; Nishiyama, Gunji & Adamatzky, 2012). Since a living system maintains a 
wholeness as a living unity and confronts perturbed surroundings, it is obliged to balance the global 
goal as a system with surrounding perturbations (Cordero, 2017; Erskine & Hermann, 2015; Gunji et 
al., 2020b). Thus, living computing is expected to solve an ill-defined problem and to obtain a quasi- 
optimal solution. Since dissipative structures such as the BZ pattern are surrounded by a specific 
boundary condition, they also contain a wholeness to some extent. In that sense, they are similar to 
quasi-biological computing. 

The third direction is the reevaluation of self-organizing criticality (SOC), in which a critical state 
characterized by a power law distribution is automatically obtained by balancing the global optimum 
with perturbation (Bak, Tang & Wiesnfeld 1987; Bak & Tang, 1989; Kauffman & Johnsen, 1991; Bak 
& Sneppen, 1993). Although there is a disadvantage for SOC that requires a fitness function (or 
potential function) that is difficult to find, the idea of criticality could contribute to the attribute of a 
quasi-optimal solution under perturbed conditions. Since it has been found that biological systems 
usually reveal the features of critical phenomena, they are expected to be systems showing SOC (e.g., 
Levi walk; Reynolds, 2018; Viswanathan et al. 1999; Bartumeus, 2007; Sims et al., 2008; Humphries 
et al., 2010; Reynolds, Schultheiss & Cheng, 2013). It leads to the idea that the results of bioinspired 
computing are destined to be quasi-optimal solutions if they show critical phenomena characterized 
by the power law. With respect to the power law distribution, the first, second and third directions of 


metaheuristics can be combined. This leads to the powerful idea that instead of implementing SOCs 


by using a specific fitness, bioinspired computing can play a more important role in generating 
criticality, which can reveal quasi-optimal solutions. 

Quasi-optimal solutions for an ill-defined problem could be generalized by universal and efficient 
computational ability under perturbed natural conditions (Gunji & Uragami, 2020). For this reason, 
universal and efficient computing devices can be made of biological materials, such as physarum gates, 
solider crab gates and some physarum computers. The next question arises (Conrad, 1983; 1985): what 
is the essential property in bioinspired or biological computing? Bioinspired computing or biological 
computing is expressed as a multi-agent system whether it is a multicellular or unicellular organism 
under indefinite conditions (Erskine & Hermann, 2015; Cordero, 2017; Erskine & Hermann, 2015; 
Gunji et al., 2020b). Therefore, these systems are obliged to be probabilistic and/or asynchronous in 
updating due to indefinite and perturbed environments. Conversely, bioinspired and biological 
computing is robust due to the probabilistic and/or asynchronous property. Thus, it is interesting to 
study the relationship among quasi-optimal solutions, critical properties, and probabilistic and/or 
asynchronous updating in a computational system (Gunji, 1990; Gunji 2014a, b; Gunji & Uragami, 
2020). This leads to the possibility of probabilistic logic gates implemented by asynchronously 
updating computations. It is important to determine the core in bioinspired computing that is robust 
under perturbed conditions. 

To answer those questions, we take a special cellular automaton called game of life (GL) (Gardner, 
1970; Berlekamp, Conway & Guy, 1982; Nordfalk & Alstrom, 1996; Rennard, 2002). Originally, 
game of life was updated in a synchronous fashion and could be used as a logic gate by using a specific 
pattern called gliders and glider guns (Berlekamp, Conway & Guy, 1982; Rennard, 2002). Since these 
patterns are generated in a synchronous fashion, their behaviors can be strictly controlled by initial 
and boundary conditions. Although the trajectories of gliders going straight can be predicted, they can 
be easily broken by small perturbations. The question of how we could construct robust logic gates 
made of gliders arises. If possible, one can implement robust probabilistic logic gates rather than 
unstable deterministic logic gates. Thus, we aim to determine the core in bioinspired computing that 


is robust under perturbed conditions 


2. Asynchronous Game of Life and its Logic gates 


2-1 Phase transition and criticality 


Game of life (GL) is a two-dimensional cellular automaton in which the state of a cell is either 0 or 


1 (Gardner, 1970; Berlekamp, Conway & Guy, 1982; Nordfalk & Alstrom, 1996; Rennard, 2002). 
The state of a cell at the (i, j) site at the tth step is represented by aj j; To compare asynchronous GL 
with synchronous GL, we introduce the virtual state of a cell at the (i, j) site at the ¢th step, which is 


represented by bi j; The transition rule of the state of a cell is defined as follows. 


In the case of aij = 0: 
If SÉ; = 3, then bf}t = 1; (1) 


Otherwisebj;* =0, 
In the case of aij =1: 
If Si; =5or 6, then bjj* = 1; (2) 


Otherwisebj;* =0, 


where 


t — wm=tin=+1 t t 
Sij F Dienst lamien aij. (3) 
The transition of synchronous GL is defined by 
t+1 t+1 
aij = bij : (4) 


In contrast, the transition of asynchronous GL is defined by using probability p with 0<p<!1 (e.g., Fatés 
2001; 2006; 2014) such that 


t+1 


aij = aij with P; (5) 


aj} = bił" with 1-p. 
It is easy to see that there are both active and inactive (stationary) phases in asynchronous GL and that 
a cell in the stationary phase is not obliged to follow transition rules (1)-(3). If p=0, then asynchronous 
GL coincides with synchronous GL. If p=1, then all cells are always at stationary phase and preserve 
given initial states over time. 

Figure 1 shows the snapshots of the pattern generated by asynchronous GL controlled by the 


probability p in equation (5), where the system size is 50 by 50. It is easy to see that under small p, the 


states converge to frozen states that oscillate with period 2 at most, while under large p, the states 


chaotically change. One can see the phase transition of the generated pattern with respect to p. 
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Figure 1. Snapshots of asynchronous GL controlled by the probability p. States 1 and 


0 are represented by black squares or blanks. 


Figure 2 shows the phase transition from ordered patterns to chaotic patterns with respect to the 
probability of the stationary phase. Patterns are generated in 500 by 500 squared cells with periodic 
boundaries. It is determined whether all cells converge to oscillatory states with a period equal to or 
smaller than 2 within 10000 steps in each trial. If true, it is counted as a frozen state. In 100 trials, the 
normalized number of frozen states is interpreted as the probability of the frozen state. Each circle 
represents the probability of the frozen state plotted against the probability of the stationary phase, p. 
To estimate the position of the critical state, the range between 0.095 and 0.155 was intensively 


investigated. The red curve represents a sigmoidal curve fitted to data, such as 


1 


v= Fae (6) 


where a = 115.038 and b = 0.1269. The critical value is approximately 0.13. Once the probability of 
the frozen state drops to 0, it is not recovered. If p is larger than 0.2, the states of asynchronous GL are 


perpetually changed in time and space. 


Probability of Frozen State 


Figure 2. The probability of the frozen state plotted against p (the probability of the 


stationary phase). The red curve represents the approximation by the sigmoidal curve. 


The next question arises, does the middle point between two phases, order (oscillation) and chaos, 
indicate the critical state. Asynchronous cellular automata showing cluster patterns or class IV patterns 
frequently reveal a power law distribution with a density of 1 plotted against time. Additionally, in 
asynchronous GL, the density against time is evaluated. 

Figure 3 shows the density of the state of 1 against time in a log-log plot. The initial condition is set 
for the density of states of 1 at 0.590. The system size is 150 by 150 cells and p=0.13. The density 
against time is averaged over 20 trials and is represented by the green curve in Figure 3. The data are 


approximated by a power-law function such as 


y = a(x —b)° (7) 


where a = 0.020 and b = 0.081, and the exponent of the power is -0.1595. 

The exponent -0.1595 could imply special meaning since this value can be found in various contexts 
(Domany & Kinzel, 1984; Hinrichsen, 2000; Fatés & Morvan, 2005). The model of directed 
percolation is defined by a squared two-dimensional system consisting of N by N cells. The state of a 
cell is either 0 or 1, where 1 means rocky material and 0 means a pore. The number of pores divided 
by N? represents the porosity. Initially, water is given in cells at the top layer. Vertically water falls 


through pores. If the water reaches at least one cell at the bottom layer, percolation succeeds. The 


normalized success rate of percolation is obtained for multiple trials, where the distribution of pores 
is randomly distributed for each trial. Directed percolation shows the phase transition between 100% 
success and 0% success of percolation. It is numerically obtained that the critical value of the porosity 
is 0.705 and that the density of water at the critical point is decayed in a power law fashion whose 


exponent is -0.1595. 
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Figure 3. The density of 1 plotted against time in a log-log scale. Data averaged over 

20 trials are represented by green curves, and the approximated line (power law) is 


represented by a purple line. 


Some rules of elementary cellular automata updated asynchronously also show the phase transition 
with respect to the probability of stationary phase (Fatès & Morvan, 2005). At the critical point, the 
density of the state of 1 decays in a power law fashion whose exponent is -0.1595. If asynchronous 
updating is defined by the order of updating, where a map of the order of cells to the order of updating 
is a bijection, and if the order of updating can influence a transition rule, many elementary cellular 
automata show cluster-like or class IV-like patterns mixing the localized periodic patterns with chaotic 
patterns. Most of them show the power law decay of the density of the state of 1 with an exponent of 
-0.1595 (Gunji 2014a, b). Therefore, this exponent, -0.1595, might show universal properties at the 


critical point in the phase transition. 


2-2 Computation by asynchronous GL 


Synchronous GL can be used as a computing resource since mesoscopic patterns called gliders can 
be used as a transmission of the value. As shown in the left diagram of Figure 4, if a special 
configuration called a glider gun is prepared, a glider that goes straight while vibrating is perpetually 
generated. Indeed, if a glider is controlled, whether it is active (i.e., gliders can be generated) or 


inactive (1.e., gliders cannot be generated), one can implement the AND gate and NOT gate. 


Va 


A=0 B=0 
| 


J 


Figure 4. Glider and AND gate made of gliders in synchronous GL. For input gliders 
A and B set at the solid squares, output A AND B (represented by A&B) is set at the 


white square. 
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Figure 5. Gliders under the perturbation of the inversion of 0 and 1, with a probability 
of 0.0001. 


The right diagram of Figure 4 shows the implementation of the AND gate by gliders. A dummy 
glider gun represented by C, hidden in the gate device, can constantly generate a glider and is prepared 
at the location in which input gliders from A and B can be collided by the dummy glider from C. If 
two gliders collide with each other, the gliders disappear. If firing a glider represents a value of 1, and 
no glider represents a value of 0, one can prepare the AAND B gate, as shown in the right diagram of 
Figure 4. This implementation reveals that if the input glider guns and the location of the output are 
stably controlled, synchronous GL has the computation ability. However, the gliders are unstable 
against perturbation. 

Figure 5 shows the behavior of gliders under the perturbation when the state of a cell is inverted 
with a probability of 0.0001 (i.e., 0 is replaced by 1 and vice versa), where there is no perturbation 
before ¢ = 200. Once a perturbation occurs, the glider never goes straight and is collapsed anywhere. 


Instead of normal gliders, broken gliders are generated and propagate radially. 
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Figure 6. Rare gliders from the initial condition where O and 1 are randomly 


distributed. Gliders are surrounded by red squares. 


Even if there is no perturbation, gliders are not stable from the initial condition, where 0 and 1 are 
randomly distributed. Gliders can be generated very rarely from the random initial condition, and they 
are collapsed by fragile broken gliders propagating irregularly. Figure 6 shows the general case soon 
after the random initial condition, in which only two gliders are generated. After this situation, local 
configurations fall into the oscillatory state called flickers. Under no perturbation, only strict and rigid 
control of gliders can implement logic gates in synchronous GL. In other words, strict and rigid 


computations such as synchronous GL cannot contribute to real-world computation, which always 


confronts perturbed conditions. Although theoretical computer scientists focus on the implementation 
of universal Turing machines or complete machines, one should think about robust but not highly 


efficient or universal computation. 
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Figure 7. Occurrence of gliders in asynchronous GL. In each diagram, gliders within 


100 steps from the initial condition overlap. 
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Figure 8. Probability of glider per cell plotted against the probability of the stationary 


phase. The peak is found close to the critical point. 


The next question arises: how can probabilistic computation be implemented against perturbed 
conditions? Here, we implement a probabilistic logic gate made of broken gliders generated in 
asynchronous GL. First, we estimate how many gliders are generated in the range of phase transitions 
of asynchronous GL. Figure 7 shows the relationship between the occurrence of gliders and the 
probability of the stationary phase in asynchronous GL. 

It is easy to see that the occurrence of gliders is increased in the range surrounding the critical value, 
p~0.13. Far from the critical value, few gliders are generated. This implies that an asynchronous GL 
with critical properties can be used as a probabilistic logic gate since many gliders that break and can 
be irregularly propagated are generated around the critical value. Figure 8 shows the normalized 
occurrence of gliders plotted against the probability of the stationary phase. The number of gliders per 
time step is normalized by the system size, which is regarded as the probability of gliders per cells. 
The curve represents the approximation by polynomial functions such as y = —3.1924x? + 


2.8931x3 — 0.9137x? + 0.1109x + 0.0034. 
2-3 Logic gate in asynchronous GL 


The probabilistic logic gate in asynchronous GL is implemented as follows. The square space 
consisting of N?cells is divided into four areas. The first area S(A) is for input value A, the second 
area S(B) is for input B, the third area S(C) is for output value, and the fourth area is the rest of the 


N?cells. These areas are defined by: 


S(A) = {6 Dli = 1,2, a2; j = 12,5}, (8) 

se) ={Gp)iaSS+1..m =o 241.9}, (9) 

sc) ={@plizt+154+2..2-1jst+154+2..2-1) co 
In the AND gate, C is indicated by A AND B (sometimes represented by A&B), and in the OR gate, C 
is indicated by A OR B (sometimes represented by A||B). In the probabilistic logic gate, the value 1 is 
defined by a density of 1 that exceeds the specific probability, P;. When a value of 1 is given for input 
X with X=A, B, each cell in S(X) has the state of 1 with the probability exceeding P4. A value of 0 for 
input X is defined by the situation in which all cells in S(X) have a state of 0. After setting input values, 
asynchronous GL with the probability of stationary phase p = 0.13 (critical value) is adapted to all 
cells in T time steps under the condition with perturbation in which state of a cell is inverted with the 


probability Pyorsg. If the density of S(C) exceeds P, at T, then output C is 1; otherwise, C is 0. Note 


that although the definition of input 0 is consistent with that of output 0, they are different from each 


other. 
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Figure 9. Typical time development of AND gate in asynchronous GL. The white 
arrow represents time development. The inner squares A, B, and A&B represent S(A), 


S(B), and S(A&B), respectively. 


Figure 9 shows the computational process of the AND gate in asynchronous GL. Here, we define T 
= 1500, N=100, P,=0.1 and pyorsg = 0.0001. This process implies that A=1 and B=1, thus A&B=1. 
The top left diagram represents an initial condition, where A=1 and B=1 are set with a density of 1 
with a probability of 0.5, which exceeds Pı. After that, asynchronous GL is adapted to all cells, with 
the boundary conditions of which are that all cells have state 0. As time proceeds, broken gliders 
generated in S(A) and S(B) are propagated to the outside of S(A) and S(B), and finally (t = 1500 = T), 
the density of 1 for cells in S(A&B) exceeds P,. This implies A&B = 1. 

Figure 10 shows typical results of the AND gate in asynchronous GL. All parameters are set to be 
the same as those in Figure 9. Since the probabilistic logic gate shows the truth table with a specific 
probability, this result is just one typical example. Each pair of patterns connected by a blank arrow 
represents a pair of initial conditions and a final pattern at t=T. This implies that if A=0 and B=0, then 
A&B = 0 (top left in Figure 10); if A=0 and B=1, then A&B =0 (top right in Figure 10); if A=1 and 
B=0, then A&B = 0 (bottom left in Figure 10); and if A=1 and B=1, then A&B =1 (bottom right in 


Figure 10). Thus, the results in Figure 10 are consistent with the truth table of the Boolean AND gate. 


The performance of this AND gate is discussed later. 

While the implementation of the OR gate is the same as that of the AND gate, if the perturbation is 
different from the condition of the AND gate, the behaviors of the logic gate are totally changed. If 
Pwoise Of the AND gate, 0.0001, is replaced by 0.001, the OR gate can be constructed. We abbreviate 
A OR B as A||B. Since broken gliders can be grown with a ten times larger perturbation (Pyoisgz = 
0.001), a pair of input 1 and input 0 can generate much more broken gliders that can reach the output 
area, S(A||B), than the situation with pyojsg = 0.0001. Thus, the probability of cells in S(A||B) with 
state | can exceed P,, which implies that the output is 1. In contrast, in the initial condition, A = 0 and 
B=0, since there is no broken glider, but instead a cell with the state of 1 surrounded by cells with the 
state of 0, such isolated cells with the state of 1 cannot be grown to the broken gliders that can 


propagate to S(A||B). That is why the result implies that output is 0. 


o 1 


Figure 10. Computational results of the probabilistic AND gate implemented by 


asynchronous GL with perturbation, Pyojsg = 0.0001, P,=0.1 and p = 0.13. 


Figure 11 shows the computational results of the OR-gate implemented by asynchronous GL with 
the probability of stationary phase p=0.13. The input and output conditions are the same as those in 
the AND gate, although only the perturbation is different, suchas Pyorsg = 0.001. As shown in Figure 
11, A=0 and B=0 lead to A||B=0; A=0 and B=1 lead to A||B=1; A=1 and B=0 lead to A||B=1; and A=1 
and B=1 lead to A||B=1. Due to a large perturbation, broken gliders can reach S(A||B). Cells with the 


state of 1 initially set as an input of 1 are first decreased and then are grown and propagated radially 


as if they had exploded. Thus, broken gliders are propagated anywhere, and the cells in state 1 could 
be grown in S(A||B). 


(0) 1 


Figure 11. Computational results of the probabilistic OR gate implemented by 


asynchronous GL with perturbation, Pyojsg = 0.001, P,=0.1 and p = 0.13. 


Now, we evaluate the performance of the AND gate and OR gate implemented by asynchronous 
GL. The conditions of the two gates are the same as those of the gates in Figures 10 and 11. For each 
input condition of a pair of A and B, it is determined whether the number of cells with the state of 1 in 
S(A&B) or S(A||B) exceeds P, (i.e., output =1) or not (output = 0) when ¢ = T. The number of trials is 
100 for each determinant, and the probability of output 1 is obtained as the number of instances of 
output | divided by the number of trials. 


Figure 12 top left, four histograms show the frequency distribution of the cover of S(A&B). They 


show histograms for the conditions of A4=0 and B=0, A=0 and B=1, A=1 and B=0, and A=1 and B=1. 


The horizontal axis represents the number of cells in S(A&B) whose state is 1 at T (1.e., the cover of 
S(A&B)), where the actual number is obtained as the number multiplied by 20 (i.e., the cover of 3 
represents 60). The vertical axis represents the frequency corresponding to the cover. The red line in 
the histogram represents P, multiplied by the number of cells in S(A&B), which implies that the cover 
on the right-hand side of the red line leads to an output of 1 and that on the left-hand side leads to an 
output of 0. Therefore, the summation of the frequencies on the right-hand side of the red line divided 


by the number of trials (100) is the probability of the output of 1. This results in the probability of the 


truth table corresponding to the AND gate, as shown in the top right of Figure 12. The probability of 
the truth table is expressed as a histogram whose horizontal axis represents a pair of inputs (00 
represents A=0 and B=0; other abbreviations of numbers are similar), and the vertical axis represents 
the probability of the output of 1 corresponding to the input pair. The histogram shows that the 
probability of the output of 1 exceeds 0.5 only for the input of 11. This result is consistent with the 


truth table of the Boolean AND gate. 


AND_00 AND_01 AND_10 AND_11 AND gate 
120 45 40 18 0.9 
m 40 | 35 16 E 0.8 
7 100 a | ; i4 a nF 
A 80 30 12 3 06 
5 25 10 2 05 
a 20 | 8 2 0.4 
2 15 2 
ù 40 15 6 = 03 
| > 2 ai 
0 0 0 Pi 0 I à 
1 3 5 7 9111315 1 3 57 9111315 1 3 57 9111315 1357911315 o 
Cover of S(A&B) Cover of S(A&B) Cover of S(A&B) Cover of S(A&B) 00 01 10 11 
OR_00 OR_01 OR_10 OR_11 OR gate 
120 18 - 16 25 12 
a 
x | 16 14 | + 
100 20 3 1 
A A 
z 80 | 12 10 i 8 0.8 
-a 10 ba 
S 60 : 8 © o 
o td 1 D 
g 40 | 6 i | | = 04 
(i i | 4 oa 
a oO 
20 a G gs 
| ; WY 7 
0 0 | 0 a 
1 3 57 9111315 1 3 57 9111315 1 3 57 9 111315 1 3 9 11 13 15 Q 
Cover of S(A| | B) Cover of S((A|| B) Cover of S((A|| 8) Cover of S((AI| B) 00 01 10 11 


Figure 12. Frequency distribution of the cover of the output area (S(A&B) or S(A||B)) 
for each input pair (left four histograms) and the probability of the output of 1 against 
the input pair (right histogram). The top corresponds to the AND gate, and the bottom 


corresponds to the OR gate implemented by asynchronous GL. 


The bottom left part of Figure 12 shows the frequency distribution of the cover of the output area, 
S(A||B), for inputs of 00, 01, 10, and 11 in the OR gate implemented by asynchronous GL. The red line 
is the same as that in the histogram for the AND gate. It is easy to see that most trials for inputs 01, 
10, and 11 exceed the red line, indicating that P, is the threshold value for the output of 1. This leads 
to the histogram for the probability of the truth table corresponding to the OR gate, as shown in the 
bottom right of Figure 12. It shows that the probability of the truth table of the OR gate is consistent 
with that of the Boolean OR gate. Finally, we can conclude that both AND and OR gates can be 


implemented by asynchronous GL at the critical point. 


3. Discussion 


In the Introduction, we discuss the metaheuristic approach based on bioinspired computing, 
unconventional computing and self-organized criticality. It is known that those approaches can obtain 
quasi-optimal solutions for an ill-defined incomplete solution and that they frequently show critical 
properties characterized by the power law distribution. How can one manipulate a search method to 
obtain a quasi-optimal solution against perturbed conditions? This question is immediately related to 
the essential core of biological and/or unconventional computing. 

Since there is no global clock in a biological system, interactions between components in a 
biological system are destined to be asynchronous. Although asynchronous computing is ubiquitous 
in natural computing systems, it has not been studied much, since its behaviors are different from 
synchronous computing, which has been extensively studied and is far from controllable. That is why 
implementing logic gates in asynchronous computing could be a touchstone to manipulating biological 
and/or unconventional computing. 

In this paper, we take two-dimensional cellular automaton, game of life (GL), and study the behavior 
of GL in asynchronous updating. Although GL in synchronous updating can generate controllable 
gliders that can be used to implement logic gates, GL in asynchronous updating cannot generate 
controllable gliders. Although gliders generated in asynchronous GL are also unstable and rapidly 
collapse into fragments, they are regenerated due to the perturbation, and broken gliders can be 
propagated. Thus, those probabilistic gliders can be used to implement probabilistic logic gates. 

First, we find that asynchronous GL shows a clear phase transition between the phase of oscillatory 
states and the phase of chaotic states. The parameter controlling phase is the probability of the 
stationary phase in asynchronous updating. In particular, we find that the critical point shows power 
law decay with an exponent of -0.1595. Since this exponent is found in various critical phenomena, it 
can be considered a universal property. 

Second, we estimate the relationship between the computability and the criticality in asynchronous 
GL. Since logic gates in synchronous GL are constructed using gliders, the computability is evaluated 
with respect to the number of gliders generated. We find that most gliders can be generated in the 
region close to the critical point in the asynchronous GL space parameterized by the probability of the 
stationary phase. This suggests that critical phenomena have autonomous computability and/or self- 
organizing reachability to quasi-optimal solutions. It is well known that most biological systems show 
a power law distribution with various properties and that living biological systems are adapted to given 
environments. Based on those results, it can be said that systems showing critical properties have the 


ability to compute fitness in environments, which can make systems themselves able to adapt to the 


environments. 

Finally, we construct the probabilistic logic gates implemented by asynchronous GL with a critical 
point. Since the generation and trajectories of gliders cannot be controlled in the logic gates, 
manipulating gliders is destined to be probabilistic. Manipulating logical values for logic gates is much 
more probabilistic than manipulating logical values in chemical computations such as BZ reaction 
computers. However, one can construct both AND gates and OR gates in asynchronous GLs. 

One of the most intriguing aspects of probabilistic logic gates implemented by asynchronous GL is 
tuning noise in manipulating logic gates. While the essential mechanism and most of the conditions 
of the OR gate are the same as those of the AND gate, only the level of perturbation can make the 
difference between the AND gate and the OR gate. In other words, increasing perturbation can replace 
the AND gate with the OR gate. This result be an essential property of living systems, since the AND 
gate outputs the value 1 only for restricted input pair 1 and 1, while the OR gate outputs the value 1 
for various input pairs, 1 and 0, 0 and 1, and 1 and 1. If the value 1 is interpreted as the active state of 
organisms, replacing the AND gate with the OR gate implies that replacing the active condition with 
the restricted conditions or the active condition with various conditions. In other words, it suggests 
that organisms utilize various resources if the environmental condition worsens (i.e., more 
perturbation). In this sense, the change from the AND gate to the OR gate by increasing perturbation 
is more biological than computational. If one can utilize the change from the AND gate to the OR gate 
corresponding to the change in the situation, a robust self-organizing computing could be implemented. 

Since the performance of the logic gate is probabilistic, distributive laws such as A&(B||C) = (A&B)|| 
(A&C) cannot always hold. It implies that the distributive law is weak. This could be related to the 


orthomodular lattice corresponding to quantum logic (Gunji et al., 2020a). 


4. Conclusion 


Derived from the idea of metaheuristic and/or self-organizing criticality, we investigate cellular 
automaton as a system carrying the essential core of biological and/or unconventional computing. 
Since the possibility of metaheuristics can be expressed as manipulating computations and the 
possibility of self-organizing criticality can be expressed as a power-law distribution, we investigate 
the behavior of an asynchronous game of life (GL) and estimate whether logic gates can be 
implemented in asynchronous GL. 

We find that asynchronous GL shows a phase transition and that the GL at the critical point could 


have the highest computability with respect to the ability to generate gliders. Finally, we implement 


probabilistic AND gates and OR gates in asynchronous GL at the critical point, which can perform as 
logic gates to some extent. In our logic gates, perturbation plays an essential role in computing. Ten 
times multiplied perturbation makes the AND gate become the OR gate under the same conditions. 
Our study can be a touchstone to spell out the novel significance of biological and/or unconventional 


computing by showing the role of noise tuned to operate logic gates. 
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